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Abstract 

A detailed study of various distinguished limits of the Green-Kubo for- 
mula for the self-diffusion coefficient is presented in this paper. First, an al- 
ternative representation of the Green-Kubo formula in terms of the solution 
of a Poisson equation is derived when the microscopic dynamics is Marko- 
vian. Then, the techniques developed in lfT5l [3ll are used to obtain a Stieltjes 
integral representation formula for the symmetric and antisymmetric parts of 
the diffusion tensor. The effect of irreversible microscopic dynamics on the 
diffusion coefficient is analyzed and various asymptotic limits of physical in- 
terest are studied. Several examples are presented that confirm the findings 
of our theory. 



1 Introduction 

The two main goals of non-equilibrium statistical mechanics are the derivation of 
macroscopic equations from microscopic dynamics and the calculation of trans- 
port coefficients El [H [3H. The starting point is a kinetic equation that governs 
the evolution of the distribution function, such as the Boltzmann, the Vlasov, the 
Lenard-Balescu or the Fokker-Planck equation. Although most kinetic equations 
involve a (quadratically) nonlinear collision operator, it is quite often the case that 
for the calculation of transport coefficients it is sufficient to consider a linearized 
collision operator and, consequently, a linearized kinetic equation. In this case, it 
is well known that the transport coefficients are related to the eigenvalues of the 
linearized collision operator |4, Ch. 13], (5J Ch. 10]. The goal of this article is 
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to present some results on the analysis of transport coefficients for a particularly 
simple class of kinetic equations describing the problem of self-diffusion. 

The distribution function f(q,p,t) of a tagged particle satisfies the kinetic 
equation 



where q, p are the position and momentum of the tagged particle and Q is a linear 
collision operator. Q is a dissipative operator which acts only on the momenta and 
with only one collision invariant, corresponding to the conservation the particle 
density. 

The macroscopic equation for this problem is simply the diffusion equation for 
the particle density p(p, t) = f f(p, q, t) dp [36] 



where the components of the diffusion tensor D are the transport coefficients that 
have to be calculated from the microscopic dynamics. 

At least two different techniques have been developed for the calculation of 
transport coefficients. The first technique is based on the analysis of the kinetic 
equation (11.11 ) and, in particular, on the expansion of the distribution function in an 
appropriate orthonormal basis, the basis consisting of the the eigenfunctions of the 
linear collision operator Q JU Ch. 13], Ch. 10]. Transport coefficients are then 
related to the eigenvalues of the collision operator. The second technique is based 
on the Green-Kubo formalism ll22l . This formalism enables us to express transport 
coefficients in terms of time integrals of appropriate autocorrelation functions. In 
particular, the diffusion coefficient is expressed in terms of the time integral of the 
velocity autocorrelation function 



The equivalence between the two approaches for the calculation of transport coef- 
ficients, the one based on the analysis of the kinetic equation and the other based 
on the Green-Kubo formalism has been studied j34l . The Green-Kubo formalism 
has been compared with other techniques based on the perturbative analysis of the 
kinetic equations, e.g. 113211311 1711 and the references therein. Many works also ex- 
ist on the rigorous justification of the validity of the Green-Kubo formula for the 
self-diffusion coefficient fISl 11011811571 151. 



-£+p-V q f = Qf, 



(1.1) 




(1.2) 




(1.3) 
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Usually the linearized collision operator is taken to be a symmetric operator 
in some appropriate Hilbert space. When the collision operator is the (adjoint of 
the) generator of a Markov process (which is the case that we will consider in 
this paper), the assumption of the symmetry of Q is equivalent to the reversibility 
of the microscopic dynamics ll33l . There are various cases, however, where the 
linearized collision operator is not symmetric. As examples we mention the lin- 
earized Vlassov-Landau operator in plasma physics Eq. 13.6.2] or the motion 
of a charged particle in a constant magnetic field undergoing collisions with the 
surrounding medium [5J Seq. 11.3]. It is one of the main objectives of this paper to 
study the effect of the antisymmetric part of the collision operator on the diffusion 
tensor. 

In most cases (i.e. for most choices of the collision kernel) it is impossible to 
obtain explicit formulas for transport coefficients. The best one can hope for is 
the derivation of estimates on transport coefficients as functions of the parameters 
of the microscopic dynamics. The derivation of such estimates is quite hard when 
using formulas of the form (11.31) . In this paper we show that the Green-Kubo for- 
malism is equivalent to a formulation based on the solution of a Poisson equation 
associated to the collision operator Q. Furthermore,we show that this formalism is 
a much more convenient starting point for rigorous and parturbative analysis of the 
diffusion tensor. The Poisson equation (cell problem) is the standard tool for cal- 
culating homogenized coefficients in the theory of homogenization for stochastic 
differential equations and partial differential equations II301L 

The problem of obtaining estimates on the diffusion coefficient has been stud- 
ied quite extensively in theory of turbulent diffusion-the motion of a particle in a 
random, divergence-free velocity field ||26l . In particular, the dependence of the 
diffusion coefficient (eddy diffusivity) on the Peclet number has been investigated. 
For this purpose, a very interesting theory has been developed by Avellaneda and 
Majda (2l|3l, see also fl"5l l6l. This theory is based on the introduction of an ap- 
propriate bounded (and sometimes compact) antisymmetric operator and it leads 
to a very systematic and rigorous perturbative analysis of the eddy diffusivity. This 
theory has been extended to time-dependent flows HI. 

In this paper we apply the Majda-Avellaneda theory to the problem of the 
derivation of rigorous estimates for the diffusion tensor of a tagged particle whose 
distribution function satisfies a kinetic equation of the form (ll.lb . We study this 
problem when the collision operator is the Fokker-Planck operator (i.e. the L 2 - 
adjoint) of an ergodic Markov process. This assumption is not very restrictive 
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when studying the problem of self-diffusion of a tagged particle since many dis- 
sipative integrodifferential operators are generators of Markov processes [19]. We 
obtain formulas for both the symmetric and the antisymmetric parts of the diffu- 
sion tensor and we use these formulas in order to study various asymptotic limits 
of physical interest. 

The rest of the paper is organized as follows. In Section |2] we obtain an al- 
ternative representation for the diffusion tensor based on the solution of a Poisson 
equation and we present two elementary examples. In Section [3] we apply the 
Majda-Avellaneda theory to the problem of self-diffusion and we study rigorously 
the weak and strong coupling limits for the diffusion tensor. Examples are pre- 
sented in Section |U Conclusions and open problems are discussed in Section [5] 



2 The Green-Kubo Formula 

In this section we show that we can rewrite the Green-Kubo formula for the dif- 
fusion coefficient in terms of the solution of an appropriate Poisson equation. We 
will consider a slight generalization of (11.11 ). namely we will consider the long-time 
dynamics of the dynamical system 

| = V W , (2.1, 

where z is an ergodic Markov process state space Z, generator C and invariant 
measure 7r(efe)Q The kinetic equation for the distribution function is 

^ + V(z) ■ V x f = C*f, (2.2) 

where C* denotes the L 2 (i?)-adjoint of the generator C The kinetic equation dl.ll ) 
is of the form (12.21 ) for V(p) = p, and where we assume that the collision operator 
(which acts only on the velocities) is the Fokker-Planck operator of an ergodic 
Markov process, which can be a diffusion process (e.g. the Ornstein-Uhlenbeck 
process in which case ( 12.21 ) becomes the Fokker-Planck equation) a jump process 
(as in the model studied in iTTll ). or a Levy process. 

Proposition 2.1. Let x(t) be the solution of (12.11 ). let z(t) be an ergodic Markov 
process with state space Z, generator C and invariant measure ir(dz) and assume 



'We remark that the process z can be x itself, or the restriction of x on the unit torus. This is the 
precisely the case in turbulent diffusion and in the Langevin equation in a periodic potential. 
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z 



that V(z) is centered with respect to p(dz), 

V(z)n(dz) = 0. 

Then the diffusion tensor (11.31 ) is given by 

D = j V(z)®<j)(z)fi(dz) (2.3) 

where (ft is the solution of the Poisson equation 

-£<p = V{z) (2.4) 

Proof. Let e by an arbitrary unit vector. We will use the notation D e = De-e, x e = 
x ■ e. The Green-Kubo formula for the diffusion coefficient along the direction e is 

r+oo 

D e = / {x e (t)x e {0))dt 



o 



(V e (z(t))V e (z(0)))dt. (2.5) 



We calculate now the correlation function in (I2.5I ). We will use the notation z = 
z(t;p) with z(0;p) = p. We have 

(V e (z(t;p))V e (z(0;p)))= [ [ V e (z)V e (p)p(z,t;p)^dp)dz, (2.6) 

J Z J z 

where p(z, t;p) is the transition probability density of the Markov process z which 
is the solution of the Fokker-Planck equation 

^ = £*p, p(z,0;p)=5(z-p). (2.7) 
We introduce the function 

V e (t,p) :=EF e (z) = J V e (z)p(z,t;p)dz 
which is the solution of the backward Kolmogorov equation 

^ = £F e , V e (0, P ) = V e (p). (2.8) 
We can write formally the solution of this equation in the form 

V e = e ct V e (p). 



We substitute this into (12.61 ) to obtain 

(V e (z(t;p))V e (z(0;p))) = / {e ct V e (p)) V e (p) fi(dp). 



z \Jo 



We use this now in the Green-Kubo formula ( 12.51 ) and, assuming that we can inter- 
change the order of integration, we calculate 

r+oo 

D e = (e ct V e (p))V e {p)n(dp)dt 
Jo 

e ct V e {p)d^j V e {p)n{dp) 
= J ((-£)- l V e {p))v e {p)n{dp) 
= Jj e V e v(dp). 

where <j) e is the solution of the Poisson equation —C(j) e = V e . In the above cal- 
culation we used the identity = J +o ° e ct ■ dt [30, Ch. 11], HU Ch. 
7]. 

□ 



From (12.31) it immediately follows that the diffusion tensor is nonnegative def- 
inite: 

D e :=e-De = [ V e <p ■ e (i(dz) = [ {-C)(t> e <t> e fi(dz) 
Jz Jz 
> 0, 

since, by definition, the collision operator is dissipative. 

When the generator £ is a symmetric operator in L 2 (Z; n(dz)), i.e. the Markov 
process z is reversible [33), the diffusion tensor is symmetric: 



Dij = / Vi(z)(j)j(z) n{dz) = I {-C)4>i{z)<j>j(z) fi(dz) 
Jz Jz 

4>i{z)(-C)4>j(z) fi(dz) = Dji. 



z 



Green-Kubo formulas for reversible Markov processes have already been studied, 
since in this case the symmetry of the generator of the Markov process implies that 
the spectral theorem for self-adjoint operators can be used ll20llT8l . Much less is 
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known about Green-Kubo formulas for non-reversible Markov process. One of the 
consequences of non-reversibility, i.e. when the generator of the Markov process 
z is not symmetric in L 2 (Z; (i{dz)), is that the diffusion tensor is not symmetric, 
unless additional symmetries are present. The symmetry properties of the diffusion 
tensor in anisotropic porous media have been studied in ||2~T1 . see also |29l . A 
general representation formula for the antisymmetric part of the diffusion tensor 
will be given in the next section. 

2.1 Elementary Examples 

The Ornstein-Uhlenbeck process. The equations of motion are 

q = p, (2.9) 
p = - 7 p + a/27/3- 1 W. (2.10) 

The equilibrium distribution of the velocity process is 

n(dp) = \ — e~^ p ' dp. 

V 27T 

The Poisson equation is 

-C4 = p, C = -jpdp + 7/3 -1 dp- 
The mean zero solution is 

<t> = -p- 

7 



The diffusion coefficient is 



D = J 4>PLi(dp) = 



7/3 

which is, of course, Einstein's formula. 

A charged particle in a constant magnetic field. We consider the motion of 
a charged particle in the presence of a constant magnetic field in the z direction, 
B = Be%, while the collisions are modeled as white noise (51 Ch. 11]. The 
equations of motion are 



% = fipxea-i/p + v^/S-VW, (2.12) 
at 
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where W denotes standard Brownian motion in R , v is the collision frequency 
and 

mc 

is the Larmor frequency of the test particle. 

The velocity is a Markov process with generator 

c = n( P2 d pl -pxdpj + vi-p-Vp + p- 1 ^)- (2.i3) 

The invariant distribution of the velocity process is the Maxwellian 

A^ViiPi 2 



Kdp) = ( e ~ m d P- 

The vector valued Poisson equation is 

-C<j> = p. 

The solution is 

v 17 17 v 1 

Pi + 9 , ^9 ^2, 9— FTjPI + 9 , ^9 ^2, -P3 



The diffusion tensor is 











u 2 +n 2 


n 


u 




u 2 +n 2 









// v 2 +iV u z +iV " \ 

p®MdV) = r l \ -jaftp ^ . (2.14) 
\ n n 1 / 



Notice that the diffusion tensor is not symmetric. This is to be expected, since the 
generator of the Markov process (12.13b is not symmetric. 



3 Stieltjes Integral Representation and Bounds on the Dif- 
fusion Tensor 

When the Markov process z is reversible, it is straightforward to obtain an inte- 
gral representation formula for the diffusion tensor using the spectral theorem for 
the self-adjoint operators |[20l . It is not possible, in general, to do the same when 
z is a nonreversible ergodic Markov process. This problem was solved by Avel- 
laneda and Majda [3 ] in the context of the theory of turbulent diffusion by intro- 
ducing an appropriate bounded, antisymmetric operator. In this section we apply 
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the Avellaneda-Majda theory in order to study the diffusion tensor (12.31 ) when z is 
an ergodic Markov process in Z. 

We will use the notation L 2 ^ := L 2 (Z; fj,(dz)). We decompose the collision 
operator C into its symmetric and antisymmetric part with respect to the L 2 ^ inner 
product: 

C = A + jS, 

where A = —A* and S = S*. The parameter 7 measures the strength of the 
symmetric part, relative to the antisymmetric part. The Poisson equation (12.4b . 
along the direction e, can be written as 

- (A + jS)^ e = V e . (3.1) 

Our goal is to study the dependence of the diffusion tensor on 7, in particular in 
the physically interesting regime 7 <C 1. 

let (•, -) M denote the inner product in L 2 ^. We introduce the family of seminorms 

2 := (/,(-S) fc A, 



Define the function spaces H k := {/ G L 2 ^ : \\f\\k < +00} and set k = 1. Then 
|| • ||i satisfies the parallelogram identity and, consequently, the completion of H 1 
with respect to the norm || • || 1 , which is denoted by H, is a Hilbert space. The inner 
product (•, ■) in % is defined through polarization and it is easy to check that, for 

f,hen, 

(f,h) = (f,(-S)h) fl . 

A careful analysis of the function space % and of its dual is presented in J24j. 

Motivated by Q, see also (6j[16l, we apply the operator (— S)^ 1 to the Poisson 
equation (13.11 ) to obtain 

(-G + 7IW = V e , (3.2) 

where we have defined the operator Q := (— S)~ 1 A and we have set V e := 
(— <S) -1 y e . This operator is antisymmetric in H: 

Lemma 3.1. The operator Q : H — > % is antisymmetric. 

Proof. We calculate 

(Gf,h) = J(-S)- 1 Af(-S)h f i(dz) = J AfhfM(dz) 

= - f (-Sr\-S)fAhn(dz) = - J f{-S)Qhii{dz) 
= ~{f,Qh). 
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□ 



We remark that, unlike the problem of turbulent diffusion [3]|6l|27]], the oper- 
ator Q is not necessarily bounded or, even more, compact. Under the assumption 
that Q is bounded as an operator from % to H. we can develop a theory similar to 
the one developed in Q. The boundedness of the operator Q needs to be checked 
for each specific example. 

Using the definitions of the space %, the operator Q and the vector V we obtain 

Dy = (3-3) 

We will use the notation || • ||% for the norm in %. It is straightforward to 
analyse the overdamped limit 7 — > +00. 

Proposition 3.2. Assume that Q : % — > T~L is a bounded operator. Then, for 7, a 
such that < 7» the diffusion coefficient admits the following asymptotic 

expansion 

D = kn 2 n + f:^\\0 k n 2 n- (3-4) 

' k=i ' 

In particular, 

lim 7D = \\V\\l. (3.5) 

Proof. We use (13.21) . the definition of the space %, and the boundendness and 
antisymmetry of the operator Q to calculate 




□ 
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From (13.51 ) we conclude that the large 7 asymptotics of the diffusion coefficient 
is universal: the scaling D e ~ - is independent of the specific properties of A, S 
or V e . This is also the case in problems where the operator Q is not bounded, such 
as the Langevin equation in a periodic potential 1 17]. 

Of course, the expansion (13.41) is of limited applicability, since it has a very 
small radius of convergence. This expansion cannot be used to study the small 
7 asymptotics of the diffusion coefficient. The analysis of this limit requires the 
study of a weakly dissipative system, since the antisymmetric part of the generator 
A represents the deterministic part of the dynamics, whereas the symmetric part 
S the noisy, dissipative dynamics. It is well known that the dynamics of such a 
system in the limit 7 — > depends crucially on the properties of the unperturbed 
deterministic system lfT3l[T4"l l9l. The properties of this system can be analyzed by 
studying the operator A. For the asymptotics of the diffusion coefficient, the null 
space of this operator has to be characterized. This fact has been recognized in the 
theory of turbulent diffusion (3]|27l[26|]. We will show that a similar theory to the 
one developed in these papers can be developed in the abstract framework adopted 
in this paper. 

Assume that Q : H -)■ U is bounded. Let TV = {/ € H : Gf = 0} denote the 
null space of Q. We have % = M ® N L . We take the projections on Af and Af 1 - 
to rewrite (13.21) as 

7<Pn = Vn, {-G + jI)cf> N ± = V N ± ■ (3.6) 

We can now write 

D e = hnt H + {ct> N ^v e NA _). 

Proposition 3.3. Assume that there exists a function p G T~L such that 

-Gp = V* ± . 

Then 

lim 7j D e = ||V£|& (3.7) 

7— s-0 

In particular, D e = 0(1/7) when = 0. 

Proof. We write 4> N ± = p + ip where tp solves the equation 

{-G + 7-0V> = -IP- 
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We use tp as a test function and use the antisymmetry of Q in % to obtain the 
estimate 

U\\h<C, 

from which we deduce that \\4> N ± \\n < C and ( 13.71 ) follows. □ 

Let Q be a bounded operator. Since it is also skew-symmetric, we can write 
G = iT where T is a self-adjoint operator in H. From the spectral theorem of 
bounded self-adjoint operators we know that there exists a one parameter family 
of projection operators -P(A) which is right-continuous, and -P(A) < P{p) when 
A < ijl and P(-oo) = 0, P(+oo) = / so that 



f(T)= [ /(A)dP(A) 



for all bounded continuous functions. Using the spectral resolution of T we can 
obtain an integral representation formula for the diffusion coefficient : 

where dfi e = (dP(X)V^ ± ,V^ ± ). We can obtain a similar formula for the anti- 
symmetric part of the diffusion tensor 

A = \{D-D T ). 
In particular, we have the following. 

Proposition 3.4. Assume that the operator Q : % — > T~L is bounded. Then the 
antisymmetric part of the diffusion tensor admits the representation 



where 
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Proof. We calculate 



Aij — „{Dij Dji) 



1 

2 

1 

2 



where we have used the notation 7£ 7 = (— zT + 7/) . From the symmetry of T 
we deduce that 1Z* = (iF + 7/)" -1 . Now we use the representation formula 



to obtain 



2 
1 
2 
1 
2 
1 
2 

1 /" \dnij{\) 



r+00 

/ e -7*( e ^_ e - 4 r^ dty | )y | 
^0 

/>+oo \ 

J e-^sm(rt)dtV[,Vlj 
e~ 7 * sin(iA) dtdfiij(X) 



2Jr 7 2 + A 2 • 

□ 

Remark 3.5. The antisymmetric part of the diffusion tensor is independent of the 
projection ofV onto the null space ofQ. 

When the operator Q : T-L — >■ % is compact we can use the spectral theorem 
for the compact, self-adjoint operator r = iQ to obtain an orthonormal basis for 
the space M^. In this case the integrals in (13.81 ) and ( 13.91 ) reduce to sums and the 
analysis of the weak noise limit 7 — > becomes rather straightforward. The weak 
noise (large Peclet number) asymptotics for the symmetric part of the diffusion 
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tensor for the advection-diffusion problem with periodic coefficients were studied 
in (6l|27l. The asymptotics of the antisymmetric part of the diffusion tensor for the 
advection-diffusion problem were studied in ||29l . 



4 Examples 

The generalized Langevin equation. The generalized Langevin equation (gLE) 
in the absence of external forces reads 

q=- I j(t-s)q(s)ds + F(t), (4.1) 
J 

where the memory kernel j(t) and noise F(t) (which is a mean zero stationary 
Gaussian process) are related through the fluctuation-dissipation theorem 

(F(t)F(s))=r 1 l(t-s). (4.2) 



We approximate the memory kernel by a sum of exponentials [23], 

N 

Under this assumption, the non-Markovian gLE (14.11) can be rewritten as a Marko- 
vian system of equations in an extended state space: 

q = p, (4.4a) 

N 

p = ^XjUj, (4.4b) 



iij = -OLjUj - Xjp + J 2f3~ l otj Wj, j = 1, . . . N. (4.4c) 

This example is of the form (12. II) with the driving Markov process being {p,u\, . . .un}. 
The generator of this process is 

c = ( E Vi) | + E ( - - x ^ + rl(Xj -L) ■ 

3=1 r j=l J J 3 

This is an ergodic Markov process with invariant measure 

pijp,u) = - e -^U+^=i^, (4.5) 
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where Z = (2ir(3 iy N+1 ^ 2 _ The symmetric and antisymmetric parts of the gen- 
erator £ in L 2 (M. N+1 ; p(p, u)dpdu) are, respectively: 



N 



d d \ 

S = J2(-a jUj — + ^ a3 — 2 ) 

j=l J 3 



and 

N N 



^=(EVi)^+E(-v 5n . 



3=1 r 3=1 

It is possible to study the spectral properties of (— S)~ 1 A. However, it is easier to 
solve the Poisson equation 

— Ccp = p 

and to calculate the diffusion coefficient. The solution of this equation is 

, Afc 1 1 

9=/ TT u k+P TT 

k— 1 Z—ik=l a k 2—ik=\ a fc 

The diffusion coefficient is 

1 



D = p4>p(p, u) dpdu = {3 1 



Z~ik=l a k 



We remark that, in the limit as N — > +oo the diffusion coefficient can become 0. 
Indeed, 



lim D 



Thus, phenomena of anomalous diffusion, in particular of subdiffusion, can appear 
in this simple model. The rigorous analysis of this problem, in the presence of 
interactions, will be presented elsewhere |[28l . 

The Generalized Ornstein-Uhlenbeck Process. We consider the following SDE 
q = p, (4.6a) 



P 



(a J - 7l)p + y / 2~/(3~ 1 W, (4.6b) 
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where q, p € Mr, J = — J , a, 7 > and W is a standard Brownian motion on 
R d . 

The presence of the antisymmetric term Jp in the equation for p implies that 
the velocity is an irreversible Markov process. The generator of the Markov pro- 
cess p is 

C={aJ- 1 I)p-V p + jp- 1 ^p- (4.7) 

— — - 1 1 2 

It is easy to check that V p • ( Jpe \p\ ) = 0. Hence, the equilibrium distribution 
of the velocity process is the same as in the reversible case: 

»p(dp) = (JLY e -§\ P \ 2 dp . 

We can decompose the generator C into its L 2 (M. d ; pp{dp)) — symmetric and anti- 
symmetric parts: 

C = aA + jS, 

where A = Jp • V p and S = — p ■ V p + /3 _1 A p . 

Using the results from [25 ] (or, equivalently, the fact that the eigenfunctions 
and eigenvalues of <S are known) it is possible to show that the operator Q = 
(— .S) -1 ^ is bounded from % := H l (M. d ; pp{dp)) and the results obtained in 
Section [3] apply. For this problem we can also obtain an explicit formula for the 
diffusion tensor. 

Proposition 4.1. The diffusion tensor is given by the formula 

D = /?- 1 (-aJ T + 7/)" 1 . (4.8) 
Proof. The Poisson equation is 

-C<j) = p, 

where the boundary condition is that <j) £ (L 2 (M. d ; [ip(dp))) d and we take the 
vector field (p to be mean zero. The solution to this equation is linear in p: 

= C P , 

for some matrix C € M. dxd to be calculated. Substituting this formula in the Pois- 
son equation we obtain (componentwise) 

/ J QkiPlCik = Pi, 
k,£ 

2 Note, however, that this operator is not compact from H to H. 
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where the notation Q = —aJ + 'yl was introduced. Notice that, since Q is positive 
definite, it is invertible. We take now the (L 2 (M. d ; fi/s(dp))) d -inner product with 
p m (denoted by (•, •)«) and use the fact that (pe,p m )i3 = /3 _1 <W to deduce 



kmPik 



l,...d. 



Or, 



Q T C 



and, consequently, C = (Q T ) 1 = (—aJ T + 7/) 1 . Furthermore, 

Dij = (<Pi,Pj)p = ( C ikPk,Pj)p 

k 



from which (14.81) follows. 



□ 



The small 7-asymptotics of D depends on the properties of the null space of 
Q := (— S)~ l A or, equivalently, A. For the problem at hand, it is sufficient to 
consider the restriction of (J\f{Q)) (or N(A)) onto linear functions in p. Con- 
sequently, in order to calculate J\f{A) we need to calculate the null space of J, 
A/"(J) = {b € M d : Jb = 0}. 

As an example, consider the case d = 3 and set 






h- > 


1 


-1 





1 


-1 


-1 






J 



In this case it is straightforward to calculate the diffusion tensor: 

-a (7 + a) —a (7 



(4.9) 



D 



1 



7 (3 a 2 + 7 2 ) 



7 2 + a 2 



a) 



a (7 — q) 7 2 + o? —a (7 + a) 



7 2 + 0? 



a (7 + a) (7 — a) 

The null space of J is one-dimensional and consists of vectors parallel to £ = 
(1, —1, 1). From the analysis presented in the previous section it is expected that 
the diffusion tensor vanishes in the limit as 7 — > along directions e±£. Indeed, 
from the above formula for the diffusion coefficient we get that (with |e| = 1) 

1 



7(3a 2 + 7 2 



■(V + !e-£|V 



(4.10) 
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Clearly, when e • £ = we have 



lim D e = 0, 

7^0 



whereas when e • £ ^ we obtain 

|e-ei 2 



e 



lim 7 D 

The diffusion coefficient, as a function of 7, and for a = 1 is plotted in Figure [Q 



5 Conclusions 

The Green-Kubo formula for the self-diffusion coefficient was studied in this pa- 
per. It was shown that the Green-Kubo formula can be rewritten in terms of the 
solution of an Poisson equation when the collision operator is linear and it is the 
generator of an ergodic Markov process. Furthermore, the effect of irreversibility 
in the microscopic dynamics on the diffusion coefficient was investigated and the 
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Majda-Avellaneda theory was used in order to study various asymptotic limits of 
the diffusion tensor. Several examples were also presented. 

There are several directions in which the work reported in this paper can be 
extended. First, a similar analysis can be applied to the linear Boltzmann equation 
(i.e. for a collision operator that has five collision invariants), in order to obtain 
alternative representation formulas for other transport coefficients, in addition to 
the self-diffusion coefficient. In this way, it should be possible to obtain rigorous 
estimates on other transport coefficients. 

Second, the effect of external forces on the scaling of transport coefficients with 
respect to the various parameters of the problem can be studied: the techniques 
presented in this paper are applicable to a kinetic equation of the form 

^ + F(q) ■ V p f + p • V q f = Qf, (5.1) 

where F(q) is an external force. 

Finally, phenomena of subdiffusion (i.e. the limit D — > 0) and superdiffusion 
(i.e. D — > oo) can also be analyzed within the framework developed in this paper. 
A simple example was given in Section 0] All these problems are currently under 
investigation. 

Acknowledgments. The author thanks P.R. Kramer for many useful discussions 
and comments and for an extremely careful reading of an earlier version of this 
paper. 
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